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EM Cygni: a study of its eclipse timings 
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EM Cygni is a Z Cam-subtype eclipsing dwarf nova. Its orbital period variations were reported in the past but the results 
were in conflict to each other while other studies allowed the possibility of no period variation. In this study we report 
accurate new times of minima of this eclipsing binary and update its O — C diagram. We also estimate the mass transfer rate 
in EM Cygni system and conclude that the mass transfer is far from the critical value. The mass transfer rate determined 
from the eclipse timings is in agreement with the spectroscopically determined value. 
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1 Introduction 

Dwarf novae are binary stars consisting of a white dwarf 
and a red dwarf in which some mass is transferred from 
the red dwarf component onto the white dwarf. Depending 
mainly on the mass transfer rate and the magnetic field of 
the white dwarf they produce outbursts caused by instabili- 
ties in the accretion disk. In the case of weak magnetic field 
of the white dwarf the main factor which governs the fre- 
quency and amplitude of the outbursts is the mass transfer 
rate. There is a critical value of the mass transfer rate and it 
was investigated by Shafter et al. (1986). According to their 
results, if the mass transfer rate is higher than this critical 
value there are no dwarf nova eruptions (we have a nova- 
like variable) and if the mass transfer rate is below it we 
have a dwarf nova. If the mass transfer is nearly equal to the 
critical value we have a Z Cam-type variable (see Shafter 
et al. 1986, 2005). (For details on dwarf novae and the pro- 
cesses in them see Warner 1995). 

In the General Catalogue of Variable Stars (Samus et 
al. 2004) sixteen eclipsing dwarf novae were listed. These 
objects are important because the period can easily and pre- 
cisely be measured due to the eclipsing nature and this al- 
lows to demonstrate so far small astrophysical effects as 
mass transfer rate. The careful analysis of the eclipse tim- 
ings can reveal small orbital period changes which cannot 
be followed by radial velocity studies. For details on the 
eclipse geometry in dwarf novae see Smak (1994). 

The subject of this study is EM Cygni which is the only 
known eclipsing Z Cam-subtype dwarf novfl Therefore the 
precise determination of its mass transfer rate is extremely 
important. Its eclipses were discovered by Mumford and 



1 The eclipsing dwarf nova V729 Sgr may also belong to the Z Cam- 
subtype, but the identification of its subtype is uncertain (Cielinski et al. 
2000). 



Krzeminski (1969) and this allows to determine its mass 
transfer rate. Pringle (1975) and Mumford (1980) stated 
that the eclipse period of EM Cyg is decreasing. In con- 
tradiction with these results, Beuermann & Pakull (1984) 
concluded that there is no evidence for period change in 
EM Cyg. Herczeg (1987) also investigated the period of 
EM Cygni and found weak evidence for period change but 
he concluded that more observations are needed. It is wor- 
thy to mention that the period decrease rate determined by 
Pringle (1975) was in the order of 10~ 10 days/cycle while 
Mumford's (1980) value is in the order of 10~ n days/cycle. 
These rather contradictory results may originate from the 
fact that the precision of minimum times obtained before 
1984 (Beuermann and Pakull 1984 estimated their accuracy 
to be 0.0008 days) could be compared to the expected width 
of the O — C in such a small period variation. If the time- 
coverage of the minimum observations would be longer, this 
effect can be avoided because the accumulated period vari- 
ations are larger in the O — C than its intrinsic scatter. 

Since 1984 more than 25000 revolutions have occurred 
and the corresponding O — C value - assuming 10~ 12 days 
decreasing in every cycle - reaches —0.0036 days. Twenty- 
one years after the latest published photoelectric minimum 
in any refereed journal we decided to observe EM Cygni 
because the accumulation of the small continuous period 
changes can be sufficient for deciding whether the EM 
Cygni's period is variable or not. This is very important 
from the point of view of understanding the mass transfer 
mechanism in cataclysmic variables. 

Our aim was twofold. First, we wanted to solve the ques- 
tion of period change of EM Cyg with new observations. 
Secondly, we wanted to calculate the mass transfer rate from 
the observed period change if such a change is present. 
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2 Observations 

EM Cyg was observed at the Konkoly Observatory on four 
nights (in Cousins I-band) and at the Baja Astronomical 
Observatory on two nights (in V-band) to obtain precise 
light curves. The observations were carried out with the lm 
RCC telescope of the Konkoly Observatory located at the 
Piszkesteto Mountain Station at 964 meters above the sea 
level and with the 50cm robotic telescope of the Baja As- 
tronomical Observatory. 

The detector of the lm RCC telescope was an 
1 340 x 1300 pixels electronically cooled (T CC d = -40°C) 
CCD camera manufactured by the Roper Scientific Ins. 
The exposure time was between 5-15 seconds depending 
on sky-conditions. The readout-time was only two seconds 
therefore the duty cycle was very favourable. All frames 
were bias-subtracted. Dark correction was not applied be- 
cause a negligible dark current was measured for this expo- 
sure time. Flat-field correction was also done using dome 
flats. The detector of the 50cm robotic telescope was an 
Apogee ALTAU16 4kx4k CCD camera with a field of view 
of 42 x 42 arcminutes. Bias, dark and flat-field corrections 
were applied. To determine the raw magnitudes of stars in 
the frames, aperture photometry was performed using the 
IRAF/DAOPHOT package (Stetson 1990). We selected a 
comparison star and a check star in the frames and they 
showed a sufficient stability in each colour (standard de- 
viations of their magnitude differences is 0.003 magnitudes 
(lm RCC) and 0.007 magnitudes (50cm robotic telescope), 
respectively). Then differential photometry was done defin- 
ing the variable's brightness as Am = m var — m comp . No 
standard transformation was applied. 

3 Analysis of the O — C diagram 
3.1 Source of data 

Times of eclipses were collected from the following 
sources. 29 minima times were taken from Mumford & 
Krzeminski (1969); two from Mumford (1974); one from 
Mumford (1975); three from Mumford (1980). One time 
of minimum was observed by Beuermann & Pakull (1984). 
These are all the published minima in any refereed journal 
and none has been published since them. 

Beuermann & Pakull (1984) determined two times of 
minima from the Figure 1 of Jameson et al. (1981). Of 
course, this method does not yield an accurate minimum 
time, but we followed them and concluded that around JD 
2444105 - when the observations of Jameson et al. (1981) 
were carried out - the O — C value was about —0.0065 days 
(against the ephemeris of Mumford & Krzeminski 1969). 
(At this time the cycle number E was close to 21390.) We 
assign a one-tenth weight to this value. 

We also found two amateur minimum obser- 
vations on the homepage of 'The Astronomer' 
(www.theastronomer.org). One of them cannot be used for 



a minimum time estimation but the other one - reproduced 
in Figure 1 - was inspected and a fitting resulted a usable 
minimum time. We denoted the reference of this minimum 
time as 'Anonymous' in Table 1. 

New and reliable minima times were determined from 
our light curves with the Kwee-van Woerden (1956) 
method which means that the descending and the ascend- 
ing branches of the minimum are reflected to a certain time 
value and it searches for the time which yields the minimum 
difference between the original and the reflected curves. 
These moments of minima are presented in Table 1 . Errors 
of the times of minima can be estimated as follows: 0.0008 
days before 1980, 0.0001-0.0004 after that. These errors are 
smaller than the observed scatter of 0.002 days seen in Fig- 
ure 4 (which shows the O — C diagram). The possible ex- 
planation of this can be due to the intrinsic variations of the 
location of the hot spot in the system. Note that the eclipses 
in EM Cyg are not total but grazing ones. Some of our min- 
imum observations are presented in Figures 2-3. 

Table 1 Times of minima of EM Cygni. References: (1): 
Mumford & Krzeminski (1969), (2): Mumford (1974), (3): 
Mumford (1975), (4): Mumford (1980), (5): Beuermann & 
Pakull (1984), (6): Anonymous (see text), (7): this paper. 
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3.2 The O-C diagram 

The Julian Heliocentric Ephemeris Date (HJED) is free 
from the varying rotation of Earth. Since we need as high 
precision as we can reach we transformed all HJD values 
of the observed minima into HJED ones. Then linear and 
quadratic fits yield the following ephemeris: 

T ocl (HJED) = 2437882.8606(3)+0.29090912(l)xE(l) 
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Fig. 4 O — C diagram of EM Cyg. The 'expected' is the 
theoretically expected O — C curve calculated from the the- 
Fig. 1 The minimum observations of EM Cyg reproduced ory (Eq. (8) of Shafter et al. 1986) assuming conservative 
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case. The 'upper limit' curve shows the maximum possible 
period variation (see text). 
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Fig. 2 The I-band light curve of EM Cyg obtained by us 
at JD 2 453 650. 
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Fig. 3 The I-band light curve of EM Cyg obtained by us 
at JD 2 453 709. 



Teci(HJED) = 2437882.8606(3) + 0.29090911(4) x E 
-2(2) • 10~ 14 x E 2 (2) 

The numbers in the paranthesis show the errors in the last 
digits. 

We also calculated the corresponding x 2 values for 
both fits and found that in case of the linear fit x 2 = 1.99 
while in case of quadratic fit it is x 2 = 1-98. So the in- 
troduction of a quadratic term does not improve signifi- 
cantly the fit. However, a simple least-squares regression 
could be very sensitive to the distribution of data points 
and of errors. Therefore we determined the value of maxi- 
mum possible quadratic term in another way. We calculated 
the O — C values within a parameter-space (in the form 
Ted = To + AE + BE 2 where E is the cycle number, and 
T = 2 437 882.8590. ..2 437 882.8622, step size: 0.0001 
days, A = 0.29090910. ..0.29090915, step size: 2.4 • 10~ 8 
days, B = -1.5 • 10~ 12 ...1.5 • 10~ 12 , step size: 10~ 15 
days) and we found that the corresponding x 2 value of 
the O - C values are minimal at T = 2 437 882.8598, 
A = 0.29090914, B = -9.8 • 10" 14 . This corresponds 
to P/P < 2.32 • 10~ 12 1/days. A Monte Carlo simula- 
tion was also used to estimate the upper limit because in 
this way we could study how the errors may yield an appar- 
ent quadratic term. At the epochs of the observed O — C 
curve, 1000 O — C values were simulated for each epoch 
with random errors. The mean of these errors was as men- 
tioned in the previous section and they were assumed to be 
normally distributed. Then we had 1000 O — C curves and 
all of them were fitted by a parabola. The average of the 
quadratic terms yielded < B >= —8.0 • 10~ 14 days. This 
is similar to the result of the previous attempt. Therefore the 
value of 2.32 • 10~ 12 1/days was accepted as the absolute 
value of the upper limit for the period decrease. 
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For comparison, we plotted the O — C diagram of EM 
Cygni calculated by the linear ephemeris together with the 
theoretically expected one in Figure 4. The maximum pos- 
sible period variation is also shown in Figure 4 which was 
calculated from the above determined upper limit of P. 

Comparing our ephemeris to the previously determined 
ones (see Section 1) one can see that our result is close to 
the one of Beuermann and Pakull (1984): we found no ev- 
idence of period variation. We cannot confirm the results 
of Pringle (1975) and Mumford (1980) who found definite 
period decrease. 

4 Discussion 

4.1 Conservative case 

Since we could only give an upper limit for the period varia- 
tion only, we could guess the maximum possible mass trans- 
fer rate in the system. Assuming a conservative case (i. e. the 
angular momentum and total mass remains constant in the 
system) the mass transfer rate can be related to the period 
change via 



P dt \q 



1 



dM 2 



(3) 



Mi + M 2 dt 

(e. g. Thomas 1977) where t refers to time, q is the mass 
ratio and M 2 is the mass of the donor star. Accepting 
q = 0.88 and M x = 1.12M , M 2 = O.99M (North et al. 
2000) and P^dP/dt = 2.32 • 10~ 12 we have dM 2 /dt = 
2 • lO~ 9 M /yr. This is an upper limit for the mass transfer 
rate. 

From the study of the far-ultraviolet spectrum of EM 
Cygni, Winter & Sion (2003) found that the mass transfer 
rate in EM Cyg is between 3.2 • lO" 11 and 1.1 • 10~ 10 so- 
lar mass per year. Our upper limit is higher than their upper 
value and is not in contradiction with their results. There- 
fore the agreement between the spectrocopically observed 
mass transfer rate and the one derived from eclipse timings 
is good. 

We compared the result to the theoretical expectations. 
Shafter et al. (1986) investigated the value of the critical 
mass transfer rate for different values of orbital periods 
and white dwarf masses. Accepting Mi = 1.12M for 
the white dwarf in EM Cyg (North et al. 2000) and tak- 
ing into account that EM Cyg has an orbital period of 6.96 
hours, Eq. (8) of Shafter et al. (1986) yields dM crlt /dt = 
1.2 • 10~ 8 M Q /yr. If the mass transfer is conservative in EM 
Cyg, this system is far from this limit. In summary, the mass 
transfer rate is lower by an order than the theoretically ex- 
pected value. 

4.2 Non-conservative case 

In EM Cyg the mass is transferred from the less massive 
component to the more massive one yielding a period in- 
crease. In non-conservative mass transfer there are several 
mechanisms which decrease the period. Therefore, in total, 



the result can be a nearly constant as well as a decreasing pe- 
riod! In this case the period variation depends on the mass 
transfer rate as well as the angular momentum (J) varia- 
tions: 



P _ J 2 + 3q Mi 3 + 2q M 2 
P~ 1~ TTq Mi ~ l + q M 2 



(4) 



(see Warner 1995). Assuming that M 2 = —M\ we have 



P = J 3(1 - q 2 ) Mi 
P J q(l + q) Mi 



(5) 



We observed nearly zero period variation so we can deter- 
mine what angular momentum loss is required to neutralize 
the period increase caused by the mass transfer. For this pur- 
pose we set P = and Mi was chosen as the critical value. 
With q = 0.88, Mi = 1.12M and Mi = 1.2 • 10- 8 M Q /yr 
we get that j are to be —4.63 • 10 _17 s _1 in order to reduce 
the period variation to zero. 

The system can loss mass via stellar wind. Assuming 
that the secondary loses 10~ 14 solar mass per year via stel- 
lar wind (which is a similar value to the one of the Sun) 
and all this mass isotropically leaves the system, this leads 
P/P = -2.4 • lO-^day" 1 (we used the formula given in 
Batten 1973, p. 95). This is too low to explain the missing 
period variation. 

The rotational period of a rotating star decreasing, if the 
star loses mass by stellar wind (Iben, Fujimoto & MacDon- 
ald 1992). Due to the spin-orbit coupling this reduces the to- 
tal angular momentum of the system. From Warner's (1995) 
Eq. (9.13b) we know that 



-1.2- 10 



34 



k 2 
OA 



^iy i2 W d y ncm (6) 



In this relation k 2 is the gyroradius of the secondary. As- 
suming k 2 = 0.1 we have j/J = -2.25 • lO" 1 ^" 1 due 
to this rotational angular momentum loss. This seems to be 
too high to explain the constancy of the period. 

The gravitational radiation is negligible in such systems. 
The magnetic braking can be a more efficient mechanism. 
This results in angular momentum loss on the scale of 

j = -2.52 • 10 34 P o 1 r g 4 (/i)dyncm (7) 

(McDermott & Tamm 1989, Warner 1995). In EM Cygni 
this yields j/J = -7.57 • lO -17 ^ 1 . This is higher than the 
required value by a factor of 2, but regarding the rather large 
uncertainties of magnetic braking theories and the fact that 
we substitute P/P = instead of its correct value (which 
is too low to determine it exactly) one can conclude that this 
itself could explain the constancy of the period. 

To hold up the mass transfer some angular momentum 
loss is required (Warner 1995). Using Warner's (1995) Eq. 
(9.16) we found that this means j/J = -1.2 • 10- 16 s _1 
which is higher by a factor of 2.6 than required. 

One can assume that more than one of these mecha- 
nisms are working in this system. But the sum of them is 
higher than the required value. If the mass transfer rate is 
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equal to the prediction of Shafter et al. (1986), then the ob- 
served constancy of the period can be explained by the mag- 
netic braking itself, and the uncertainty of the estimations 
allows that we can regard the effects of mass transfer and 
the magnetic braking to be equal. Also, magnetic braking 
can be the process which drives the mass-transfer (because 
its magnitude is in the required range). But in this case Win- 
ter & Sion (2003) have been measured too low mass transfer 
rate. 

The mass transfer rate given by Winter & Sion would 
mean an angular momentum loss rate of J/J = —2.7 • 
10 _19 s _1 if the period variation is zero. This is close to 
the range of the angular momentum loss caused by the as- 
sumed value of the secondary's stellar wind rate (note that 
this real or accidental agreement does not mean that this 
should be due to the angular momentum loss caused by stel- 
lar wind) but it does not fit the rate of other kind of possible 
angular momentum loss. With their mass transfer rate the 
magnetic braking and/or the rotational braking would dom- 
inate the right side hand of Eqs. (4-5) and we would observe 
a large period decrease. One can estimate the period varia- 
tion rate from their M value. Substituting their mass transfer 
rate into Eq. (5) and assuming that there is an angular mo- 
mentum loss from magnetic braking (Eq. 7) we would have 
P/P = — 2.0 • 10 _11 days~ . From the available minima 
observations this figure should be excluded. 

5 Conclusion 

We observed six new times of minima of EM Cygni. We 
updated its O — C diagram and found that the period is con- 
stant. We found that the upper limit for the period decrease 
is ||| < 2.3- 10~ 12 1/days. 

If the mass transfer is conservative then its rate is far 
from the theoretically predicted value (Sect. 4. 1) but it is in 
agreement of the spectroscopically determined mass trans- 
fer rate of Winter and Sion (2003). (Note that their method 
was indirect and its reliability depends on e.g. the accu- 
racy of distance therefore this agreement between the two 
independent determinations give a hint only that the mass 
transfer is not conservative.) In conservative case the spec- 
troscopically observed mass transfer rate and the one given 
by present eclipse timing analysis is lower than the theoret- 
ically predicted value. 

If we assume that the mass transfer occurs in a non- 
conservative mode and the rate of mass transfer equals to 
the theoretical predictions of Shafter et al. (1986), we found 
that the period increase caused by the mass transfer from 
the secondary to the primary is reduced by the angular mo- 
mentum loss due to magnetic braking to zero (Sect. 4.2). 
This explains why we did not observe period variation. It is 
worthy to mention that in this case the spectroscopic mea- 
surement of the mass transfer rate by Winter & Sion (2003) 
is in contradiction. 

Therefore their spectroscopic measurement would ex- 
clude the scenario we investigated in Sect. 4.2. If their mea- 



surement is right then the mass transfer cannot occur in non- 
conservative case and this would mean that the accretion 
disk theory model of dwarf novae may be challenged by the 
currently available observations. 

Since EM Cyg is the only known eclipsing Z Cam- 
subtype dwarf nova its further minima observations are im- 
portant to refine our results and determine more precisely 
the difference between the observed and expected mass 
transfer rates. 
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